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Abstract 



. We study a new approach to determine the asymptotic behaviour of quantum many-particle 



systems near coalescence points of particles which interact via singular Coulomb potentials. This 
problem is of fundamental interest in electronic structure theory in order to establish accurate 
and efficient models for numerical simulations. Within our approach, coalescence points of par- 
ticles are treated as embedded geometric singularities in the configuration space of electrons. 
Based on a general singular pseudo-differential calculus, we provide a recursive scheme for the 
calculation of the parametrix and corresponding Green operator of a nonrelativistic Hamiltonian. 
^ ' In our singular calculus, the Green operator encodes all the asymptotic information of the eigen- 

functions. Explicit calculations and an asymptotic representation for the Green operator of the 
hydrogen atom and isoelectronic ions are presented. 

g ; 1 Introduction 

O . 

The asymptotic behaviour of eigenfunctions of nonrelativistic quantum mechanical Hamiltonians 
near singular points of the potential is of great physical interest and paramount importance for 
the numerical solution of these eigenvalue problems. This is especially true in the field of elec- 
tronic structure calculations where electrons and nuclei interact via singular Coulomb potentials 
with each other. We present a new approach to extract asymptotic information by an explicit re- 
cursive calculation of the parametrix of the Hamiltonian and corresponding Green operator. Such 
an approach, however, requires the full machinery of pseudo-differential calculus on manifolds with 
singularities. The Coulomb potentials in the Hamiltonian induce a hierarchy of embedded singu- 
larities into the configuration space of the electrons, i.e., M. 3N for a system of N electrons within 
the Born-Oppenheimer approximation. This means that we have to consider M. 3N as a stratified 
space with conical, edge and higher order corner singularities. The underlying general calculus of 
pseudo-differential operator algebras and corresponding function spaces has been developed by one 
of the authors in great generality, cf . the monographs [TJ [2] . It has been already applied to various 
types of boundary value problems with applications in material science, i.e., crack theory, cf. the 
joint monograph [3] by two of the authors. Applications to quantum theory are natural from a 
conceptual point of view, however, they face considerable technical difficulties with regard to the 
tremendous complexity of this calculus. Especially for higher order corner singularities it is still 
under development. It is therefore the purpose of the present work to provide an explicit calculation 
along the line of the general approach, which illustrates the basic idea and its potentiality. 

The paper is organized as follows: first we briefly summarize previous work on the electronic 
Schrodinger equation, focusing on the asymptotic behaviour near singularities of the potential. 
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In the subsequent section, we present a general outline of the basic ideas of the calculus thereby 
avoiding most of the technicalities. In order to make the paper reasonably self-contained, we present 
basic elements of the calculus in Appendix |Bj The main part of the paper is devoted to explicit 
calculations for the hydrogen atom and isoelectronic ions, where we discuss the recursive construction 
of the parametrix which enables us to extract the asymptotic information to any order. 

1.1 Previous work on the asymptotics of many-electron systems 

The rigorous mathematical analysis of the asymptotic behaviour of eigenfunctions of nonrelativistic 
Hamiltonians near singularities of the potential was initiated by Kato [JJ. Presently a considerably 
more detailed picture is available by the work of M. and T. Hoffmann-Ostenhof and coworkers 
[51 EE El El EH]- Their work provides some leading order terms of asymptotic expansions and 
presents a source of inspiration for numerical methods [IT], which try to incorporate such terms 
explicitly into the basis set. For certain purposes, however, like nonlinear approximation theory 
|12j . a more stringent control of the singular behaviour is required. In this particular case, these 
results only provide some guideline to conjecture the more general asymptotic behaviour. 

There are only a few explicitly solvable electronic Coulomb systems mentioned in the literature. 
The best known example is the hydrogen atom dicussed below. For systems with more than one 
electron, it is only for the helium atom where an explicit solution has been conjectured by Fock 
|13j . Despite important rigorous results in favour of Fock's expansion by Morgan [14j . a complete 
proof is still missing. However, separating the leading order terms of the Fock expansion enabled 
improved regularity results for the remaining part of general many-electron wavefunctions, cf. [9]. 
Unfortunately this approach is presently limited to leading order terms because spin-degrees of 
freedom are not taken into account. 

The conical singularity that occurs when an electron approaches a nucleus, which is considered 
to be fixed within the Born-Oppenheimer approximation, is the simplest case in the hierarchy of 
singularities discussed before. It is well-known that the eigenfunctions of the hydrogen atom possess 
a simple Taylor asymptotic with respect to the electron-nuclear distance. Such an asymptotic 
behaviour has been recently established for the nonlinear Hartree-Fock model [15] using tools from 
a singular pseudo-differential operator calculus. This result provides, e.g., a solid ground for the 
application of nonlinear approximation theory to singular particle models |16| . 

1.2 Fuchs type differential operators and their parametrices 

Due to space restrictions, we cannot provide a detailed account of the operator calculus on manifolds 
with conical singularities. For the convenience of the reader we have summarized some basic features 
of the calculus in Appendix [Bj It is the purpose of this section to give a brief outline of the main 
ideas underlying the present work. 

The cone algebra over an infinite stretched cone X A = R + x X for a closed compact C°° manifold 
X may be motivated by the task to establish a (pseudo-differential) operator calculus that contains 
all Fuchs type cone degenerate differential operators of the form 



for cij 6 C°°(lR+, DifP _:, (X)) , together with the parametrices of elliptic elements. Here DiS u (X) 
denotes the space of differential operators on X of order v with smooth coefficients in local coor- 
dinates. In such a context it is natural to admit /x E N rather than /j, = 2 which is the case for 
our applications, since compositions and parametrices give rise to other orders. Under some natural 
growth conditions on the coefficients a,j(r) for r — >• oo the operator A induces continuous operators 
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between weighted Sobolev spaces 

A : JC S '^(X A ) — > IC s - 2 ^- 2 (X A ), 

for all s, 7 G K. A brief discussion of these spaces can be found in Appendix IB. 11 For a detailed 
account, we refer to the monographs [IJ [2]. 

Our approach to characterise asymptotics of eigenfunctions of cone degenerate operators (II. ID 
is based on a singular calculus of reasonable generality which means that it can be applied to 
a large class of problems where the present one is only a specific example. Hereby it is a basic 
task to understand the structure of parametrices. In general, asymptotics, together with weight 
information and smoothness of solutions, may be expected from a parametrix P of the respective 
elliptic operator A. The left product PA is equal to the identity plus a smoothing operator G, where 
the latter contains all the information provided that P is known in detail close to the singularity. 
Within the present work, we restrict to the case of conical singularities with smooth base manifolds. 
More general edge and corner singularities will be discussed in a subsequent publication. 

In order to deal with parametrices it is necessary to consider a wider class of pseudo-differential 
operators. The corresponding symbols, e.g., for Fuchs type differential operators (jl.ip . are (up to 
the factor r~^) locally of the form 

PA(r,x,p,g) = p(r,x,rp,£), 

where p(r,x,p,^) is a classical pseudo-differential symbol in the covariables (p, £), smoothly de- 
pending on (r, x) up to r = 0. The Fuchs type symbol pa(t,x, £) is degenerate for r — > with 
respect to the covariable p. Globally along the base X of a cone one considers operator functions 
p(r, p) = p(r, rp) where p(r, p) is a parameter-dependent family of classical pseudo-differential opera- 
tors on X with parameter p € R, smoothly depending on r € R+, i.e., p(r, p) G C°°(R+, L^(X; R p )). 

Elliptic Fuchs type differential operators (jl.ip can be alternatively expressed as Mellin pseudo- 
differential operators along M + with holomorphic Mellin symbols h(r, z) taking values in (classical 
pseudo-) differential operators on X. This admits the construction of a parametrix with symbol 
M -1 )(r, z) representing the Mellin-Leibniz inverse of h(r, z) with meromorphic Mellin symbol which 
has poles in the complex z-plane. The poles and their multiplicities just determine the asymptotics 
of eigenfunctions near r = 0. 



2 The hydrogen atom revisited 



In nonrelativistic quantum theory, the hydrogen atom represents the most prominent textbook exam- 
ple for a quantum mechanical Hamiltonian where all eigenfunctions belonging to discrete eigenvalues 
are explicitly known. The stationary Schrodinger equation for the hydrogen atom and isoelectronic 
ions has the form 

(H-E)u = 0, (2.1) 

with Hamiltonian operator 



H 



(2.2) 



2 \x\ 

in atomic units, where A is the Laplacian in R 3 9 x and Z the charge of the nucleus. In polar 
coordinates, i.e., (S 2 ) A := R + x S 2 = Kj \ {0}, equation (J27TJ) becomes 

2 



Au 



1 



d_ 

Or 



d_ 

Or 



+ ^A 5 2 + rZ + r 2 E 



u 



0. 



(2.3) 



It can be seen that the electron-nuclear Coulomb potential can be incorporated into the Fuchs type 
differential operator acting on (S 2 ) A with embedded conical singularity at the origin. 
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Let cj,u),u) be cut-off functions such that ojCj = oj, ojCj = Cj. Throughout this paper a cut-off 
function is any real- valued oj € Cg°(R+) such that oj(r) = 1 close to r = 0. It is easy to show that 
A can be decomposed into two parts in the following way 

A = ujACj + (1 - oj)A(1 - Cj). 

Alternatively the first term can be expressed as Mellin pseudo-differential operator with holomorphic 
operator valued symbol 

h(r,w) = \{w 2 -w + A 52 ) + rZ + r 2 E, 
and A in the second term can be pushed forward to the differential operator E — H, i.e., 

A = ojr~ 2 o^ l (h)Cb + (1 - u){E - H)(l - Cj). 

In the sequel we investigate A as an operator in the cone algebra, corresponding to the stretched cone 
(5 2 ) A . (We refer, for example, to [1, Chapter 8] for details on such operators and corresponding 
symbols; a brief outline is given in Appendix IB. li ) 

As a crucial point first we consider the bijectivity of the principal conormal symbol of the Mellin 
operator 

a M (A)(w) = \(w 2 -w + A 52 ) , 

as a parameter dependent second order differential operator between standard Sobolev spaces on 
S 2 . 

Lemma 1. The principal conormal symbol defines isomorphisms 

a M (A)(w) : H S (S 2 ) -> H S - 2 {S 2 ) (2.4) 
for any s£l provided that w £ Z. 

Proof. Prom the ellipticity of a^(A) follows the Fredholm property of (|2.4p for any s € R. Therefore 
it is sufficient to consider the (of s independent) kernels of (|2,4p and its adjoint. Let us take 
s = 2. The eigenfunctions of the operator A52 : H 2 {S 2 ) — > L 2 (S 2 ) are the spherical harmonics 
\Xim '■ I € No,—/ < m < 1} corresponding to the 2^ + 1 times degenerate eigenvalues — 1(1 + 1). 
These functions form a complete orthogonal basis in L 2 (S 2 ) and belong to H oc (S 2 ). A simple 
calculation shows that the kernel of o~m(A) has no nontrivial element for w £ Z. The adjoint 
operator a M (A)* : L 2 (S 2 ) -> H~ 2 (S 2 ) satisfies (u,a M (A)*v) = (a M (A)u,v) for all u G H 2 (S 2 ). If 
v belongs to the kernel of o~m(A)* , we get (o~m(A)u, v) = for any u € H 2 (S 2 ); in particular, 

(a M (A)Y lm , v) = i (w 2 - w - 1(1 + 1)) (Y lm , v) = 0. 

It follows that if w £ Z we have {Y[ m , v) = for all l,m and therefore v = 0. □ 

Let us now check the ellipticity conditions for A (as an element in C 2 ((S 2 ) A ,g),g = (7,7 — 2, O) 
with = (—oo,0]) in the sense of [lj Definition 6, Section 8.2.5]. The first condition is obviously 
satisfied. The second condition follows for any 7 ^ Z + i by Lemma 1. It only remains to check the 
exit behaviour of A. The corresponding exit symbols 

a e (A) = -\ (t 2 + i 2 + £f) + E and a^ e (A) = -\ (£? + f 2 + , 

have constant coefficients and satisfy for E < (which is the case in our application) the required 
exit conditions, i.e., 

a e (A) ^ for all f e M 3 and a^ e (A) + for all ^Gl 3 \{0}. 

From the ellipticity of A it follows the existence of a parametrix and corresponding Green operator. 
The explicit construction of the parametrix requires the inverse of the principal conormal symbol 
&m(A) as a meromorphic operator function with poles at the non-bijectivity points wq G Z. 
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Lemma 2. The inverse of the principal conormal symbol o~m(A) 1 { w ) is a meromorphic operator 
function with simple poles at wo G Z; more precisely, 

l 

Ea M (A)' 1 (w) = -2{w + ly 1 (;Y lm }Yi m for w < and I = -w 

m=—l 
I 

Ho- M (^) _1 (w) =2(w-l- l)" 1 (^Yi m )Yi m for w > and I = w - 1 

m=— £ 

where Yi m := ^pp^ Y\ m and S denotes the principal part of the Laurent expansion of ctm(^1)~ 1 (w) ^ n 
a neighbourhood of the pole w$ . 

Proof. The lemma is an immediate consequence of Theorem 7.1 in [17]. In order to apply the 
theorem it is only necessary to calculate at w = wq the maximum of the null multiplicity of 

(w 2 - w + A 52 ) fi m (w, 0, 4>), 

for arbitrary w holomorphic functions with fi m (wo,9,4>) = Yi m (6,(p). Suppose the first derivative 
with respect to w vanishes at w = wq, i.e., 

(2w - 1) fim(w , 9, 4>) + [wq -w + A S 2) fi m (w , 9, 4>) = 0, 

which can be written as 

(tug - w + A s2 ) f( m (w , 9,4) = - (2w - 1) Y lm (9, 4>). 

This equation has no solution for f[ m (wo,6,4>) since the right side belongs to the kernel of the 
self-adjoint operator on the left side. Therefore all null multiplicities are of first order. □ 

According to the spectral theorem for operators with purely pointwise spectrum (cf. [H]), we 
can represent the inverse of the principal conormal symbol by the absolutely convergent sum 

X p I 

o'm{A)~ 1 (w) =2V- 1 with Pi := V (;Y lm )Y lm . (2.5) 

ur — w — 1(1 + 1 

1=0 \ < J m =-l 

An alternative simple proof of Lemma [2] can be given using the explicit representation (|2.5p . 
2.1 Construction of the parametrix and Green operator 

The operator A belongs to the cone algebra on (S* 2 ) A and defines bounded operators between 
weighted Sobolev spaces 

A : lC s 'mS 2 ) A ) — > JC S - 2 '^~\(S 2 ) A ). 

In the following, we require the corresponding weighted Sobolev spaces with asymptotics K, S q ((S 2 ) a ) . 
These spaces can be considered as direct sums 

/C^((SY) = ^ + /C^((Sy), (2.6) 

of the flattened weighted Sobolev spaces 

JCl?((S 2 ) A ) := p| IC s ^- e ((S 2 ) A ), 

e>0 
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with = 0], — oo < < 0, and asymptotic spaces 

L j k=0 > 

The asymptotic space Sq is characterized by a sequence qj G C which is taken from a strip of the 
complex plane, i.e., 

Qj € : | - 7 + < &s < | - 7 j > 

where the width and location of this strip are determined by its weight data (7, G) with 6 = 0] 
and —00 < $ < 0. Each substrip of finite width contains only a finite number of qj. Furthermore, 
the coefficients Cjk belong to finite dimensional subspaces Lj C C°°(S 2 ). The asymptotics of Sq is 
therefore completely characterized by the asymptotic type Q := {(pj, rrij, Lj)}j e ^ + . In the following, 
we employ the asymptotic subspaces 

^((SY ) := [a G /C^((S 2 ) A ) : (1 - u)u G 5(R, C°°(S 2 ))| R+x52 } , (2.7) 

with Schwartz class type of behaviour for exit r — > 00. The spaces /Cq 7 ((5 2 ) a ) and <Sq((S' 2 ) a ) are 
Frechet spaces equipped with natural semi-norms according to the decomposition f|2.6[) ; we refer to 
[U [2] for further datails. 

As an elliptic element of the operator class C 2 ((S 2 ) A , g), g = (7,7 — 2,0) with 7 £ Z + \ 
and = (—oo,0], A has a parametrix P in the cone algebra which belongs to C~ 2 ((S 2 ) A ,g), 
g = (7 — 2, 7, 0). It can be written in the general form 

P = Jr 2 op]7 3 (h^ (r, w))u/+(l-u/)P(l- u') , 

where a/, a)', a)' are cut-off functions satisfying u/a;' = a;', w'a)' = Co' and P represents a standard 
pseudo-differential operator of order —2 on R 3 . By definition, the parametrix satisfies the equation 

PA = 1 + G withGGC G ((S 2 ) A ,0,), ^ = (7,7,0), 

where the Green operator G maps /C S ' 7 ((S' 2 ) A ) into Sq for some discrete asymptotics Q. Application 
of the parametrix to (|2.3p yields u = —Gu, i.e., the asymptotic behaviour of the eigenfunctions 
belonging to the eigenvalue E are completely characterized by the corresponding Green operator. 

In order to proceed with our discussion, we need an educated guess for the parameter 7. Taking 
7 to large, it might happen that a solution of Au = does not exist in /C s ' 7 and a too small 
value for 7 may introduce additional solutions which are too singular to be meaningful from a 
physical viewpoint. In general it should be not too difficult to identify a priori an appropriate 
interval for 7 using standard regularity results. In particular, for the hydrogen atom explicitly 
known eigenfunctions belong to /C s ' 7 with 7 < |. Therefore we assume ^ < 7 < | in the following. 

After we have established the notion of a parametrix and Green operator in the pseudo-differential 
algebra, it is possible to state our main theorem. 

Theorem 1. The Hamiltonian of the hydrogen atom (Z = 1) and for isoelectronic ions (|2.2p belongs 
to the cone algebra, i.e., it belongs to the operator class C 2 ((S 2 ) A ,g), g = (7,7 — 2, 0), and is elliptic 
in the sense of [i, Definition 6, Section 8.2.5] for 7 ^ Z + i. There exist a corresponding parametrix 
P G C- 2 ((5 2 ) A , ff - 1 ), g 1 = (7 - 2,7,0) and Green operator G G C G ((S 2 ) A , 9l ), 9l = (7,7,6). 
The Green operator for | < 7 < | has a leading order asymptotic expansion of the form 

G(u) = u' (1 - Zr + (Z 2 - E)±r 2 + ■ ■ ■) P Q (u) 

+oj'{Zr-±(Zr) 2 + ±(Z 2 -2E)Zr 3 + ---)P 1 Q 1 (u)+ ••• (2.8) 
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with coefficients 



Qo{u) 
Qi{u) 



2 7 -l 



(M(oj" - 1) op^ 1 (h ) (urn)) (0) + (Mu"(rZ + r 2 E)u) (0) 

(M(w" - l)o P X7 1 (/io + Zr)(wu))(-1) + (Mu"r 2 Eu)(-l) 
(M(w" - 1) op^ 1 (h ) (urn)) (0) + (Mw"(rZ + r 2 E)u) (0) 



given by values of Mellin transformations M (cf. Appendix \B.1\) at non-bijectivity points of the 
principal conormal symbol ctm(A), where u £ /C S,7 ((5 2 ) A ), uj" := cj'w, and Pi, I = 0,1,2,..., 
denotes projection operators on subspaces which belong to eigenvalues —1(1 + 1) of the Laplace- 
Beltrami operator on S 2 . 



Remark 1. Using the explicit expression for eigenvalues of the hydrogen atom and isoelectronic 



ions, i.e., E n 



111 



G n (u) 



, n € N, we obtain 



J (1 - Zr + (1 - ^z)\(Zr) 2 + ■■■) PoQoi 



u 



W (Zr - \{Zr) 2 + i(l - ^){Zrf + ■ ■ ■) P 1 Q 1 (u) + ■ ■ ■ 

which can be easily verified for asymptotic expansions of the analytic expressions of the eigenfunc- 
tions for ground and excited states with angular momentum quantum number I = and 1, cf. the 
monograph [19 1 . 

In order to calculate the asymptotics of the Green operator it is sufficient to focus on a neigh- 
bourhood of the embedded conical singularity at r = 0. Using the general commutation relation 



r n op p M (g(r, w)) r~ n = op^ +n (T n g(r, w)) 

with shift operator T n acting on the Mellin symbol by T n g(r, w) = g(r, w + n), we get 

PA = u'r 2 op]^ 3 (h^~^(r, w))tJj'u}r~ 2 opj^ 1 (h(r, w))u) + ■ ■ ■ 
= w' op]^ 1 (T 2 /i (-1) (r, w))Q'u op]^ 1 (h(r, w))ui H 



(2.9) 



(2.10) 



where the omitted terms do not contribute to the asymptotics of the Green operator. For the follow- 
ing calculations it is convenient to change our notation, we define ho(w) := om{A)(w), h\(w) := Z 
and h2(w) := E. Herewith the holomorphic symbol becomes 

h(r, w) = ho(w) + rh\(w) + r 2 h2(w). 

Similarly, we perform a Taylor series expansion of the Mellin symbol of the parametrix 



N 



h(-V{r,w) = ^2r i h[- 1 \w) + R N (r,w). 



i=0 
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Inserting both expansions into (I2.10P yields 



PA 



CO 



' °Pm 1 (j2 rWht^iw) + T 2 R N (r,w)\ Co'uj op]7 1 (h (w) + rhy (w) + r 2 h 2 {w))Co + 
n=o ' 

N 

' r* op]^ 1 (T 2 ^ -1 ^ (w))u>'u oplJ 1 (h (w) + rh^w) + r 2 /i 2 H)^ + ■ ■ ■ 
i=0 

AT 

' £ op]," 1 (T 2 ^- 1 ) (u,))o/ w op]; 1 (ftoM)fi 



JV 

' Y, r l+1 opU 2 (r 1 ^ HKu, op],' 1 (fo («,)) 



i=0 
JV 



8=0 



The last two lines contain products of Mellin operators with different weights. Adding and sub- 
tracting terms with Mellin operators of the same weight yields 



N 



PA 



2=0 

N 

W £ r i+1 op],- 1 (r 1 ^^^))^ 0PI7 1 (/nH) 

i=0 
JV 

W J3 r l+2 OPX7 1 (w))u/u OPX7 1 (^H)o; 



i=0 



N 



+ 53r i+1 W 'opl7 2 (T 1 / li ( - 1) H)^-a; / opl7 1 (r 1 / l (- 1) H)^ w op^ 1 (h t (w)) u> 



i=0 
JV 



8=0 



+ ^ r*+ 2 a/ op]7 3 (fcH> (w))u/ - J op 1 ' 1 (w))u/ u> OP17 1 (M™))" 



7-l^(-l)/ 



where the operators in the last two lines will be identified as a part Gj of the Green operator G. In 
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the next step we identify the remaining part of the Green operator G 

N 



PA = W '^r J opX7 1 (T 2 ^- 1) H)op]7 1 (/ l0 H)cD 

N 

W £ op]; 1 (T 1 ^- 1 ) («,)) op]," 1 (w)) 



Co' 



i=0 
N 



W £ r*+ 2 op]," 1 («,)) op],' 1 (fo (w))w 



i=0 



W ^ r* op],- 1 (T 2 ^- 1 ) («,)) (u/ w - 1) op^ 1 (ho{w))Q 

N 

W £ r l+1 op]," 1 (T 1 /^ («,)) (£/ w - l) op^ 1 (h, («,)) 



Co' 



> =: G// 



i=0 



W ^r^ 2 opX7 1 (^ 1) H) (tf w - 1) qp^faM) 



Co' 



i=0 



+ • 



The individual operators in the first three lines and the omitted terms are classified according to 
their asymptotic behaviour for r — >■ 0. Therefore, the term corresponding to r° is equal to the 
identity operator and all other terms corresponding to r l , i = 1, 2, . . ., must be zero. This provides a 
recursive scheme to calculate the meromorphic operator valued symbols h\ ^ (u;) of the parametrix 
explicitly For r° we get 

i = op^r^-V)) op^^H) = o P X7 1 (T 2 4- 1) (^)/ l0 (^)) 

and therefore 

(T 2 4 _1) H)^oH = 1 — >■ h ( - 1 \w) = T- 2 h 1 (w) = h 1 (w-2). (2.11) 
The condition for j = 1 is 

= op^T^-V)) op^M™)) + op^T 1 /**" 1 ^)) opX7 1 (/ ll (u,)), 



or 



o = (t^V^M™) + (V/^V))/^), 

which gives us 

h[- 1] H = -[T- 1 hi- 1 \w))(T' 2 h 1 (w)){T- 2 h \w)] 
= -Zh - 1 (w-3)h Q - 1 (w-2). 

Similarly for j = 2 

o = o P ]7 1 (T 2 4- 1) H)op]7 1 (/ l0 H)+opX7 1 (T 1 / l S- 1) H)op]7 1 (/ ll M) 
+ °PM 1 ( /i o _1) H) op^^^H), 



(2.12) 
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or 



0= (r 2 h { 2 1) ( w j\h (w)+ (r l h[ 1 \w]\h l {w) + h { - 1) (w)h 2 (w), 

which together with (f2TT|) and (j27T2j) yields 



M = - (zT^h^ (w) + ET- 2 h { ~ l) (wj) T~ 2 /tQ 1 



= Z 2 /io 1 (u)-4)/io 1 (u--3)/i ( 7 1 ('u;-2) -^/i " 1 (u;-4)/io 1 (u>-2). (2.13) 

Obviously, the recursive scheme can be easily continued up to arbitrary order. For our present work 
we consider only terms up to j = 2 which is sufficient to demonstrate asymptotic properties of the 
Green operator. Using the spectral decomposition (|2.5p . we finally obtain 



Pi 



/r } H=2V ^_ ' (2.14) 



/ q {w + i_ 2){w _i_ 3 y 
^ P 



h[ 1 \w) = -azY j ( —j „ v —j J — t^tt — r^-y (2-15) 

^— ' [w + I — o)(w + I — 2)(w — I — 3)(u> — 1 — 4:) 



1=0 

4 _1) H = E f 7 ; If ; ^ " 4^ 

2 v ' \ (hi 4- Z — 3 Win — / — 4^ 



z=o 



(to + Z - 3) (to - Z - 4) j (w + Z - 4) (to + Z - 2)(w - Z - 3) (to - Z - 5) ' 

(2.16) 



Because of PjPj' = for Z 7^ Z', the poles of ^(u>) and ^(w) remain simple. 

Next, we consider Gi and G//-type Green operators in more detail. In the following let us set 
Yp := {w G C : $iw = (3} for any (3 G R. The operator contains operators of the general form 



u/op]^ (h(w))u)' — u/op^ (h{w))£o' 
J op^ 1 {Kw))r^<uj' - Jr<^ op^ 1 (T-^H)^ 



r ~7+7_ 



According to [21 Proposition 2.3.69], the term in square brackets belongs to Cg{{S 2 ) a , g) with 
g = (7,7, (—oo,0]). By repeating the calculations in the proof of the proposition mentioned before, 
in our case 7 > 7, the whole expression can be transformed into 

oj' op^ 1 (h(w))v — cj'r 7- ^ op^ M 1 (T~ 1+ ^h(w)) r~ 1+ ^v 

1 



J I r~ w h(w)(M^ 1 v)(w)dw, (2.17) 

2ttz J r -i ur 3 _ 

with u := w'u and u G /C S -T((S 2 ) A ). Since (M^_i«) (w) belongs to A$ - 7 < < § - 7, C°°(S 2 )) , 
we can apply Cauchy's theorem to replace the integral contour by a closed smooth curve Sj 7 j C 
{w G C : I — 7 < Ku> < I — 7} which encloses all poles in the strip and compute the integral using 
Cauchy's formula. 

The G//-type Green operators can be computed in a similar manner. The operators are of 
general form 

w ' °Pm 1 {9(w)) (u" - l) op]J 1 (h(w))ui, 

with oj" := Cj'u and g(w) G M^(S 2 ), h(w) G M V R {S 2 ) (cf. Appendix for some n,v G R and 
discrete asymptotics ii, Q. Because of u/u;' = uJ and wo) = a; we can assume w.l.o.g. that suppa/ n 
supp(l — lo") = 0. It follows from [21 Lemma 2.3.73] that this operator belongs to Cg{{S 2 ) a , g) with 
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g = (7, 7, (—00, 0]). Following the proof in [2] we choose a cut-off function ui with (1 — u/')(l — ui) = 
(1 — to") and observe that 

(l-w)qp^ 1 (/»(«;)) w : )C S ^((S 2 ) A ) ^ u" K s ~ v ^ +N {{S 2 ) A ) + (1 - u")H s - v ^ {{S 2 ) A ) 

is continuous for every TV € N. Furthermore 

u/opl^H) (a/' - 1) : J'K S ~^ +N {(S 2 ) A ) + (1 - a/Oft 8 ""' 7 ((S 2 ) A ) — ► S 7 fe ((S 2 ) A ) 

is continuous, where the asymptotic type P^ has weight data (7, (— k, 0]) with € N sufficiently 
large for iV. It follows from these considerations that it is the mapping 

u/opl^HKu/'-l) :H s -^((S 2 r) ^Sl((S 2 ) A ), 

which determines the asymptotics of the Green operator. For u € T-L s ~ u, ' y r (S 2 ) A ) and every 7 > 
7 it follows that (M(u" - l)u)(w) belongs to A(% - 7 < < | - 7,C°°(5' 2 )). Because of 
•H^ t/ ' 7 ((5' 2 ) A ) C W~ v ^ {(S 2 ) A ) , we can rearrange the operator into 



w '°Pa/ [g(w))(u>" -1) =u' op], f (flr(iy)) - op]./ (g(w)) (u" - 1) + J op^ (w" - 1), 

where the last term maps 

ft*-"^ ((£2)A) into /c~.7((52)A) and therefore does not contribute to the 
asymptotics for sufficiently large 7. The first term in square brackets can be computed like in the 
case of the Gi Green operator, i.e., 



oplf 1 (ff W) - oplf 1 (g(w)) (u" - 1) (it) 



,7-1 



—j\ r~ w g{w)(M^ l {oj" -l)u){w)dw. (2.18) 

2m Jr 3 ur- 1 



3 ui o 
2"^ |-7 



Proof of Theorem [TJ 

Proof. The explicit expressions for the asymptotic behaviour of the Green operator G can be easily 
calculated from (i2~Tlt (|2"7T51) . (pOol) using Eqs. (pTTjl and ([215]) . The details of the calculation are 
given in Appendix [A] Collecting individual terms order by order with respect to powers of r and 
the angular momentum projection operators Pi, I = 0,1,2..., we obtain the asymptotic formula 

(EH). □ 



Appendices 



A Calculation of Green operators 

In this appendix, we characterize Gi and G/j-type Green operators. 

A.l G/-type Green operator 

We begin with the Gi operator 



N 



Gl = Z^r l+l 



i=0 



N 



+2 



8=0 



/ 7— 3/i(— 1)\ // / 7— 1)\ // 

w °Pm W ) w - w °Pm C«i ) w 
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where we have computed contributions up to N = 2. 

Using (|2.14p and (|2.17p . we obtain for i = with v := d'u 

" ■ , qp3T a (T 1 4" 1) ) - a/qplr 1 ^^)] W 

= Zra/ 



Zr 



27ri ./r^ 1 ur 



3 Ui 5 „, 



and 



r^T 1 /^ x) («;)(M 7 _2t;)(«;)d«> 
= -2Za/fb(M«)(l), 

^«p2r 8 (4" 1) )- w '^pjr 1 (4" 1) )l («) 

r- w h { ~ 1} (to) (M 7 _ 3 t;) («;)dti; 
£?a/ [2fb(Mu)(2) + |rJ\(M«)(l)] . 



PrV — 



./r^ 1 ur 



1-7 4-T 



For i = 1, we get with (|2.15p 



Zr 2 



>o V \f(T^)-Jo V l~\T^)\{v) 



Zr 



2, ./ 1 / ^-WrplA-l) 



2vri ./r- 1 ur 



r ~w T i h (-i) ^ ( Mj _ 2 v) (w)dw 
Z 2 d [-2rP (Mv){l) + |rPi(Mu)(l)] , 



§-7 3-T 



and 

Er 3 



CO 



E r 3 d— 

2?" ./r- 1 ur 



r -, "4 1) (^)(M 7 _ 3 v)(w)du; 



|-7 4-T 



= -EZd [-2rP [Mv)(2) + §rPi(Mu)(2) - ir 2 Pi (Md) (1) + ±r 2 P 2 (Mi;)(l)] 
Finally, for i = 2 we get with (|2.16p 



Zr a 



= Zr 3 u 



„3, ./ 1 



2vri ./r- 1 ur 



-Zd 



r- w T l h ( 2 1) (w)(M 1 _ 2 v){w)dw 

| r 2 ( Z 2 _ ^ p Q (1) _ 1 r 2 Z 2 pi ^ Mv j ^ + l_ r 2 ( Z 2 + 2 ^ p 2 ( Mv j (1) 



a Ul 5 
|-7 S-T 



and 
Pr 4 



doplfih^-doplJ 1 ^') (v) 



Pr 4 w 



4,7 1 



2-7TZ /r- 1 



r"^4 1) (w)(M 7 _ 3 v)(w)d-u; 



r 3 ur 7 



2 r 2 / 'r 7 2 



Z 2 - P) P (Md) (2) - ir 2 Z 2 P x (Afu) (2) + ^r 2 (Z 2 + 2E) P 2 (Mv) (2) 



+ ^r 3 (Z 2 - 2E)P 1 (Mv){\) - ±r 3 Z 2 P 2 {Mv)(l) + J,r 3 (±Z 2 + P)P 3 (Mt;)(l) 
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A. 2 G/7-type Green operator 



According to our discussion in Section 12.11 the G//-type Green operator is given by 



A? 



G u 



i=0 

N 

+^E 



°V" M \T 2 ^ 1) ) ~ 0V^(T^) (u" - 1) o P r (M* 



.7-1, 



i=0 
N 

i=0 



[UJ — 1 UJ 



1)' 



where the dots indicate a remainder which does not contribute to the asymptotics. The parameter 
7 > 7 is chosen sufficiently large such that all poles belong to the intervall [| — 7, | — 7]. 

Using (|2.14p and (|2.18p . we obtain for i = with v := uju 



UJ 



oV 1 m\T 2 ^)-o^ m \t^- 1) ) {J'-l)o^-\h Q ){v) 



UJ 



2ni 



r 3 ur" 1 

5-T |-7 



v E^TT r ' p ' ( M 7-i(^" - iWST 1 (*>)(«)) H)> 

2=0 



Zruj' 



Zrw'-^- / r- w T 1 /t^ 1) (u;)(M 7 _i(a;" - 1)(«)) 



Zw 'E ^7TT r ^ ( M 7-i(^ - 1)(«)) H + 1), 



1=1 



Er 2 uj' 



17 2 / 1 

2?ri 



r-™4 -1) H ( M 7-i ( w " - !)(«)) ( w ) dw 



r 3 ur~ 

3-7 |-7 



- Eu} ' E 2[Ti r ' P; ( M ^- l(u;// " H + 2) - 



Z=2 
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For i = 1, we get with (j2.15[) 

op^T 2 /^) -opl,- 1 ^ 2 /^)] {u"-l)o V ^\h ){v) 
, 1 



2-ni 



1=0 



r 3 urT 1 



(Z + 1)(2Z + 1) 



ZrV 



- Zw ' E pTi) r ^ ( M 7-i K - 1) opIZ 1 M («)) H + 1), 
"opXr'Cr 1 ^) - qpJF 1 ^" 1 ')] - i)(f) 

= ZrV— /" r-T^" 1 '^) (M,_i(w" - 1)(«)) (tfl)d«) 

27TZ J r , 



! ur- 1 

S-7 |-7 



^(Z + l)(2Z + l) 



r l+1 P l {M^ 1 {u;"-l)(v)) (-1 + 1) 



zV E uwzr/ p i ( M 7-i(^' - 1)(«)) H + 2). 



PrV 



1(21 + 1) 

°pxr 1 c4" 1) )-opJr 1 ('4" 1) )] v-m 



/iW-^ / r- w h { f 1] {w) (Af 7 _i(w" - 1)(«)) (w)dti; 



8-t §"7 



^(Z + 1)(2Z + 1) 



r ,+1 P,(M rl ( W "-l)(«)) H + 2) 



1=3 



Finally for i = 2, we get with (|2.16p 



rV 



OPI7 1 (T 2 h^) - oplJ 1 (T 2 ^) (a," - 1) opXZ 1 (fc>) (v) 



r 2 oo 



2iri 



[ r- w T 2 h { 2 1] (w) (M 7 _i(o/' - 1) op^ 1 (h ) (vj) (w)dw 

Jr 3 ur; 1 v 7 



/r |-^§-7 

z 2 



(2Z + l)(2Z + 3) 



r'+ 2 P (^.jK - 1) op^ 1 (Zi ) («)) (-Z) 



^Z(Z + l)(2Z + l) 
Z 2 \ 2 



«=2 



r' +1 P ; (M 7 _xK - l)^ 1 (/*,)(«)) (-Z + 1) 
(2Z - 1)(2Z +"l) r ^ K-iK " 1)opX7 1 (Z 10 )(^)) (-Z + 2), 
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o^-\T^) - oJ-\T^)] (J' - l)(v) 

= Zr 3 uj'— f r-T^^H - 1)(«)) (u»)dw 

L • - 2 N 2 



(2Z + l)(2Z + 3) 



r^P^M^" -!)(*;)) (-Z + 1) 



1=3 



z 2 



E 



(21 - l)(2l + 



— r P; (M 7 _i(oj" - l)(v)) (-1 + 3), 



Er A J 



o^-\h^)-o^-\h^)] (u"-l)(v) 

= Er 4 uj'^- [ r- w 4 _1) (^) (M 7 _i(w" - l)(v)) 



3 Ul , 

j-7 f-7 



(2/ + l)(2/ + 3) 



r ,+2 fl (M 7 _!( W " -!)(«)) H + 2) 



+ EZ2 "'E i(i + 1) ( 2i + 1) ^ - DM) H + 3) 



B Tools for operators on a cone 

We complete here some details on the cone algebra that belong to the technical background of this 
paper. 

B.l Mellin pseudo-differential operators on a cone 

Before we start our discussion of Mellin pseudo-differential operators, let us briefly recall the def- 
inition of /C s ' 7 (X A )-spaces for the case of X = S n , the unit sphere in R 1+n . By the identification 
between M. 1+n \ {0} and X A via polar coordinates x — > (r, x) we get 

1C S ~<(X A ) = ujH sr/ (X A ) + (1 - u)H s (R 1+n ). 

Here n s ^(X A ) = rm s >°(X A ), and rL s '°(X A ) for s £ N is defined to be the set of all u(r,x) € 
r- n / 2 L 2 (R + x X) such that (rd r ) j Du € r - n / 2 L 2 (R + x X) for all D £ DiS s ~ j (X), < j < s. The 
definition for s € R in general follows by duality and complex interpolation. 

In the definition of pseudo-differential operators on manifolds with conical singularity it is natural 
to use the Mellin transform in the direction of the cone variable r. Let us set Tp := {w G C : $tw = (3} 
for any real (3. We systematically employ the Mellin transform 

f'°° dr 
(Mu)(w) = / r w u(r) — , 
Jo r 
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which is for u G Cq°(M.+) an entire function in w G C, where Mu\r^ G Sftp) f° r every /3 G R, 
uniformly in compact /3-intervals. M 7 u := M-uIri is called the weighted Mellin transform of 

weight 7 G R. We also admit functions/distributions in r G R+ with values in a Frechet vector 
space, e.g., Cq°(R + , C°°(X)). Then M 7 , first applied on elements with compact support in r G R + 
extends to more general spaces. For instance, when we define H s (Fp x X) to be the Fourier transform 
along of the space H s (Fp x X) (the usual cylindrical Sobolev space of smoothness s G R) then 

M 7 _n : ojK s ^(X a ) # s (r«+i x X) 

'2 2 ' 

is continuous for every s G R. 

After these preparatory remarks, we can now introduce the concept of Mellin pseudo-differential 
operators on X A . Let V\ (X, R z ) denote the space of classical parameter-dependent pseudo-differential 
operators of order p G R on X, with parameter A G R . This is a Frechet space in a natural way. 
When R is identified with for some j3 we also write L"~JX, Fp). For any operator- valued amplitude 
function f(w) G L^(X, Fn+i T ), we obtain a Mellin pseudo-differential operator 

op], 2 (/)wu := M-_ x f [f{w)M 1 _n U ju(w)). 
Multiplied by another cut-off function Co from the left this gives us a continuous operator 

wop]; ! (/)w : K S ^(X A ) -> 1C S -^(X A ) 

for every s G R. It is known that such a continuity also holds for / = f(r, w) G C°° (R+, L^(X, Fn+i )) . 
An example is 

i=o J 

compare with the formula The corresponding polynomial in w G C gives rise to a holomorphic 

function with values in DifP(X), smoothly depending on r G R + . Let Mq(X) for any |i£K denote 
the subspace of all h £ „4(C,L^(X)) such that /i|r« 6 L^ t (X, Fp) for every /3 G R, uniformly in 
compact /3-intervals. The space Mq(X) is Frechet in a natural way, and it makes sense to talk about 
C°°(R+,M£(X)). In particular, we have Y.^=o a j( r ) wj G M e>PO- 

Let us denote by cr^(A) the homogeneous principal symbol of A of oder /x as an invariantly 
defined function on T*X A \ 0. Then in the variables (r, x) and covariables (p, £), we also have, what 
we call the reduced principal symbol, 

a^{A)(r, x, p, £) = r'V,/, (A) (r, x, r~ x p, £). 

By ellipticity of ^4 with respect to <r^ (briefly, cr^-ellipticity), we understand the condition a^{A) ^ 
on T*X A \ and a^(A)(r, x,p,£) / for (p, £) ^ up to r = 0. In this case, we always know that 
the principal conormal symbol 

a M (,4)H := J^Oj(0)t^ : # S (X) -> F S -^(X) (B.l) 
i=o 

is a family of Fredholm operators, and that there is a sequence (pj)jeZ C C such that — > oo 
as |j| — ?> oo, and (|B.lj) bijective for all u; G C \ (j>j)j £ z- Now A is said to be elliptic with respect to 
a weight 7 G R if (j)j)jez H Fn+i_ 7 = 0. Then (jB.ip is a family of isomorphisms for all w G f n +i 
Furthermore, if in addition A is exit-elliptic (cf. [2]; some details are also outlined in [15]) which is 
a condition for r — > 00 then 

A : /C S ' 7 (X A ) /C S -^^(X A ) (B.2) 
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is a Predholm operator for all s € R. Actually, it can be proved that the ellipticity conditions are 
even necessary for the Fredholm property of (|B.2p , 

It is a basic task (not only for the specific application in the present paper) to assess asymptotics 
of solutions to an equation Au = f for r — > when / G /C S ~^' 7 ~^(X A ) has an asymptotic expansion 
of the kind 

rrij 

f(r, x)~^2J2 d jk (x)r-** log k r (B.3) 

j k=0 

for r — > 0, with coefficients djk £ C°°(X). In the case of our application we have / = which is, of 
course, a special case. It turns out that any solution u(r, x) has an analogous asymptotic expansion 
as (|B.3|h where the points p € (pj)j^z from the system above with < — 7 contribute further 
exponents to the resulting asymptotic expansion. The generalities may be found in [2]. Here we 
sketch a few details, mainly to explain notation that played a role before. 

B.2 The structure of parametrices 

Asymptotics of solutions may be regarded as an aspect of elliptic regularity, here in subspaces 
/Cp 7 (X A ) C )C S,J (X A ) of functions u(r,x) with asymptotics of type P, i.e., 

rrij 

u(r, x) ~ c jk( x ) r ~ Pj l°£ fc r 

3 k=0 

as r — > 0. The asymptotic type P = mj)}j £ ^ C C x No is a sequence of pairs where $lpj — > —00 
as j —> 00, and Cjk € C°°(X). In the notion of asymptotics, we may also include the coefficients 
and control the finite dimensional subspaces of C°°(X). For simplicity, we drop this aspect in the 
general description. 

The weighted Mellin transform M 7 _« , n = dimX, transforms LoJCp" 1 (X A ) for any cut-off function 
w to a space of functions that are meromorphic in w € C with poles at pj of multiplicity rrij + 1. This 
picture fits to Mellin operators with meromorphic symbols. Given a sequence Q := {(qj, nj^jez C 
C x No with \ ffiqj\ — > 00 as \j\ — > 00, by Mq°°(X) we denote a suitable space of functions f(w) € 
^4(C\-7rc<5, L~°°(X)) for ircQ := {qj ■ j € Z} which are meromorphic with poles at qj of multiplicity 
rij + 1. We call Q a Mellin asymptotic type. Similarly as for functions we can also control the 
coefficients in the principal parts of the Laurent expansions to be finite rank smoothing operators 
on X; however, again for simplicity, we ignore this aspect in our discussion. The space Mq°°{X) is 
defined by the condition that for every / £ Mq°°{X) and vrcQ-excision function \ 

for every /3 6 R, uniformly in compact /^-intervals. Moreover, we set 

M*(X) := M£(X) + M Q °°(X). (B.4) 

A basic observation is that when A is elliptic with respect to oy,, there is a Mellin asymptotic type 
Q such that 

a M {A)- l (w)eM Q ^{X). 

This is just what we exploit in our paper. Another aspect of our method is to use the existence of 
a parametrix P of A which is, roughly speaking, of the form 

P = opj^ ^ 2 (h) + remainder 

for a function h(r, z) € C°°(lR + , Mq^(X)) with a suitable Mellin asymptotic type Q. The remainder 
mainly consists of a so-called Green operator. Such an operator is defined by its mapping properties, 
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namely, to map /C s_/i ' 7_/t (X A ) to a space of the type Sj,(X A ) for a suitable asymptotic type P (and 
a similar condition for the formal adjoint). In the more "orthodox" cone algebra that may be found, 
e.g., in [2], one refers to a decomposition of h similarly as (|B.4h . according to 

C°°(R + ,M Q M (X)) =C 00 (R + ,M > M (X)) +C 00 (R + ,M Q 00 (X)), 

and localises the Mellin operators near by cut-off functions on both sides. What we do here is to find 
the coefficients of the Taylor expansion of h(r, w) for r — > which belong to Mq fJ '(X) and to compose 
operators localized at 0, cf. (|2.1U|) . where we employ commutation properties between operators with 
meromorphic Mellin symbols and powers of r, cf. (|2.9p . which leaves Green remainders. The latter 
ones are just explicitly computed in our example. 
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